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B : ABSTRACT 

The lifetime of whispering gaUery modes in a dielectric cavity with a metallic in- 
clusion is shown to fluctuate by orders of magnitude when size and location of 
the inclusion are varied. We ascribe these fluctuations to tunneling transitions be- 
^ ^ tween resonances quantized in different regions of phase space. This interpretation 

' is confirmed by a comparison of the classical phase space structure with the Husimi 

■ distribution of the resonant modes. A model Hamiltonian is introduced that de- 

o ■ . 

\^ , scribes the phenomenon and shows that it can be expected in a more general class 

i of systems. 
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In recent years a lot of experimental effort pl has been devoted to dielectric 



O 

^ . 

Cj . microresonators. In these resonators long-lived whispering gallery (WG) modes are 
created by total internal reflection of light circulating just inside the surface of the 
I dielectric. WG modes can have Q values exceeding p 10^ and they are of great 
' interest for applications such as microlasers and optical interconnects. For optical 
^ cavities that are deformed from rotationally symmetric shape it was demonstrated |@ 

that the emission pattern of WG modes can be understood in a ray-optics picture. 
The ray dynamics typically exhibits both stable and chaotic trajectories, and it was 
argued that the presence of chaos should lead to a broadening of the WG modes with 
increasing deformation. 

In this paper we study effects in dielectric cavities that are beyond a ray-optics 
model. We show that tunneling transitions between classically disconnected regions 
in phase space can lead to fluctuations in the lifetime of WG modes by several orders 
of magnitude. This is demonstrated numerically by solving the wave equation for a 
dielectric that has a permeable coating and the shape of the annular billiard. The 
fluctuations arise from avoided crossings of WG modes with much broader resonances 
of the dielectric. The coating serves to make this effect more pronounced by raising 
all lifetimes above those obtained from total internal reflection alone. The Husimi 
projections of the relevant modes reveal that the broad resonances and the WG modes 



are localized in chaotic and regular regions of phase space, respectively. This moti- 
vates us to model the dielectric in terms of a tunneling Hamiltonian where states 
quantized in different regions of phase space are connected by small tunneling matrix 
elements. We derive the quantum scattering matrix associated with the tunneling 
Hamiltonian and obtain the widths of the quasibound states. In agreement with our 
numerical observation we find strong fluctuations in the lifetime of WG modes close 
to avoided crossings with modes localized in a different region of phase space. 

The idea that disconneted regions in classical phase space can give rise to 
quantum tunneling has previously been introduced by Davis and Heller [^], who 
dubbed this notion dynamical tunneling. Dynamical tunneling has been observed 
B S 0) m ill the spectrum of closed quantum systems whose classical counterparts 
exhibit both regular and chaotic motion. In such systems dynamical tunneling leads 
to statistical fluctuations in the splittings of nearly degenerate energy levels. While 
the resulting level splittings are typically so small that their experimental detection 
seems unlikely, we believe that the tunneling effects giving rise to fluctuations in 
resonance lifetimes of open systems are amenable to experimental verification. 

We numerically study a dielectric cavity with the shape of the annular billiard 
0, ||. The system [Fig. 1 a)] consists of a cylindrical dielectric with radius R and 
index of refraction n > 1 that is suspended in air. Embedded in this cylinder and 
off-center by an amount 5 is a totally reflecting (metallic) rod of radius a < R — 6. 
The surface of the dielectric is covered by a penetrable metallic film. The dielectric 
function of the system is written as e = e*^*^^"^ + e^^™ where e^*^^*^ = inside the 
dielectric and e^"^"^ = 1 outside. The metallic film is taken as a (5-function layer, 
gfiim _ —fjfi}Ji^[jr — R^ with ?7 > 0, corresponding to a purely imaginary refractive 
index. This neglects any phase shifts in the film (caused by a real part of the index), 
as well as the frequency-dependence of the absorption. Absorption by the inclusion 
is neglected because we expect the dominant loss to be caused by the outer interface. 

We solve for the quasibound states of the annular billiard using wave function 
matching. The electric field \I' when polarized parallel to the cylinder axis satisfies 
the scalar wave equation 

V2^ + e(r)P^ = 0, (1) 

and is continuous at r = i?. The 5-function in e creates a jump in the derivative, 
d'^/drlR = ?7n^ k^R'^{R). We search for a solution with no incoming waves, 

^<=E^M*M' ^>=T.B,H^^\kr) COS ^^<t>, (2) 
^+ = Hf\nkr) cos/i0+^ Si^^H^^\nkr) cosi/0, (3) 

where \l/< and \&> denote the electric field for r < i? and r > R, respectively. We 
use polar coordinates r, with origin in the center of the dielectric, and H^^-^^ are 
the Hankel functions of first and second kind. The scattering matrix S"^^-* of the inner 
annulus is known explicitly p[. Matching ^ and d'^/dr at the surface r = R, one 



0.05 0.1 0.15 

(5/R 



a) Cross seection of cylindrical resonator with ofF-ccntcr metallic inclusion (inner disk), b) Width of 
the WG mode (2,30) vs. eccentricity 5/R. Inset: Positions of WG mode (2,30) (solid line) and the 
closest broad resonance (dashed). 



obtains a set of simultaneous equations for the Fourier coefficients A^, i?^, which can 
only be solved at discrete complex values of the wavenumber k. To each solution we 
assign the radial and angular momentum quantum numbers (m, m) of the state in the 
concentric billiard {5 = 0) from which it evolved. 

The c/osec? annular billiard {t] oo) exhibits mixed classical dynamics 0. Let 
X be the angle of incidence measured from the surface normal. Then WG trajectories 
are those which never hit the inner cylinder, i.e. their impact parameter is R \ sin x\ > 
a+6. Their motion is regular, with a conserved x- The phase space for R \ sin x\ < d+S 
contains both regular and chaotic regions if 5 > 0. For chaotic rays, x fluctuates 
between reflections and can become very small. In the open system, this means 
that modes supported by the chaotic domain have widths that are several orders of 
magnitude broader than the widths of WG modes, because the outer coating is far 
more penetrable close to normal incidence than at the grazing angles allowed for WG 
rays. 

Fig. 1 b) shows the width of the even WG mode (2,30) corresponding to 
sinx = 0.66 (at 5 = 0) as a function of 6/R. We chose n = 2, t] = 1 and kept 
a + 6 = 0.65 fixed. This ensures that WG modes with | sinx| > {a + S)/R remain 
classically undisturbed. The width is seen to fiuctuate by orders of magnitude. This 
is caused by avoided crossings with broader resonances, as we show in the inset for 
the large peak at 6/R ~ 0.139 0. Here, the dashed line corresponds to a broad 
resonance of width —lm{kR) ~ 10^'^. 

We now want to associate these resonances with different regions of the clas- 
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The projection psos of the Husimi function onto the Poincarc surface of section, superimposed with 
classical trajectories to indicate the extend of regular and chaotic regions. The uncertainty in (j) is 
0.5 radians; 6/R = 0.135 in a) and b), with a) showing the phase space density of the WG mode 
and b) that of the broad resonance. In c), psos is shown for the WG state at S/R = 0.139. 



sical phase space spanned by r, and their conjugate momenta {Pr,P(j))- This is 
achieved using the Husimi function, i.e. the overlap of \l/< with a minimum-uncertainty 
wavepacket $ centered at f, 0, pr,P(f,- The resulting phase space density Ph{^, 4>, Pr, P<i>) 
characterizes \E'<. To obtain a classical phase-space portrait, we employ a Poincare 
surface of section with and sin x as coordinates. By semiclassical arguments, sin x 
can be directly related to p^. One can then project pn onto the surface of section by 
integrating out Pr, choosing the spread of $ around f to be zero and setting f = R. 

The resulting density distribution psos is shown in Fig. 2 for the two reso- 
nances singled out in the inset to Fig. 1 b). Away from the avoided crossing, the 
WG mode is localized in the regular regime at sinx 0.66, whereas the broad res- 
onance is supported in the chaotic domain. However, close to the anticrossing at 
~ 0.139, Psos [Fig. 2 c)] reflects a strong coupling between the regular and 
chaotic parts of phase space due to interference of the two resonances. In contrast to 
Figs. 2 a) and b), this is a classically forbidden phase space distribution made possible 



only by dynamical tunneling. It causes the lifetime of the WG mode to approach that 
of the chaotic state. We find fluctuations similar to Fig. 1 b) for all WG modes with 
I sinxl > {a + 5)1 R. 

Our analytical approach is motivated by the observation that the modes of 
the cavity have overlap with distinct regions of the classical phase space. Exploiting 
the analogy between the two-dimensional wave equation and the two-dimensional 
Schrodinger equation we model the system by the following Hamiltonian 

H = Ho + Ht. (4) 

Here, Hq describes states quantized in two classically disconnected regions in phase 
space as well as the channel region when these subsystems are totally disconnected 
from each other. The couplings between the subregions are described by Hrp- Written 
in a basis where the various subsystems are diagonal, Hq has the form 

i/o = E ^^lh + E E.clc^ + E / dEEdlEd^E. (5) 

Here, the creation (annihilation) operators qj and {qi and c^) arise from the quan- 
tization of the two respective regions of phase space. The corresponding quantum 
states will be referred to as regular and chaotic states below. A generalization of Hq 
to include more regions in phase space is straightforward. There is a continuum of 
channel states with corresponding operators d'l^, d^^, where a = 1, 2, . . . , M denotes 
the channels. The tunneling Hamiltonian Ht is given by 

Ht= (E / dEUa^{E)diEq^ + H.c) + 

\ ai / 

+ (e + if.c. j + (^E / dEW,^{E)d\,,c^ + H.c^ . (6) 

It couples both the regular states and the chaotic states directly to the continuum. 
Moreover, there is a coupling between the regular and the chaotic states via the 
matrix elements V^^. 

The Hamiltonian previously has been investigated in various limiting cases. 
The case V^^ = Uai = together with the assumption that the matrix elements Wafi be 
random variables defines a model for chaotic scattering which has been the subject of 



intensive research in the past decade 0, ^ . The case of vanishing channel couplings 



Uai = Wa/j. = has been studied in recent work [4-6] on the splitting distribution in 
closed quantum systems. Here, we shall investigate the effects of tunneling transitions 
between the regular and the chaotic region {Vi^ ^ 0) when these regions are weakly 
coupled to the channels. The channel couplings are assumed so weak that the typical 
widths of both the regular and the chaotic resonances are smaller than their respective 
energy spacing (the regime of isolated resonances). 



We calculate the scattering matrix Sab{E) = 5ab — '^'^iTab at energy E from 
the Lippmann-Schwinger equation 



T = Ht + Ht^^—-—-T (7) 
E - Ho + 17] 

for the transition operator T (t] is positive infinitesimal). This equation is solved 
by the explicit summation of the associated Born series |]12|. Here it is sufficient to 
discuss the solution in the case that the direct coupling between regular and channel 
states is much weaker than the coupling arising from indirect processes involving 
chaotic states. One finds that the S matrix has the form Sab = ^ab — 27riT^ — 2mT^ 
where comprises all scattering processes involving only chaotic states while all 
remaining processes are included in T^. The poles of are located at complex 
energies Ei = ti — {i/2)Ti^ where is the position and Fj the width of the i-th 
resonance. Assuming that this resonance is isolated from all other resonances, the 
result is 
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where the matrix elements are evaluated at energy Cj. This resembles the expression 
for the transition rate in second-order perturbation theory, but in contrast to the 
perturbative expression the denominator in Eq. has a finite imaginary part. In 
accordance with our numerical results for the lifetime of WG modes, Eq. (|^) shows 
a large increase of the width whenever a regular WG mode (with index i) crosses a 
chaotic mode (with index fi). In this way, multi-step tunneling processes involving 
chaotic states cause strong fiuctuations in the resonance widths. 



It is known |13] that the interference of resonances may strongly affect their 
positions and widths. To study this effect very close to an avoided crossing we consider 
the nonhermitian Hamiltonian 




as a simple model of two resonances that are coupled to single open channel. If the 
separation \Si — E^\ is much larger than \V\, there are two eigenstates with width 
close to and T, respectively. For small values of \£i — E^\, one finds an avoided 
crossing of the resonance positions (the real parts of the eigenvalues). At the same 
time, both widths approach the value r/2. We expect the sharp resonance typically 
to aquire a width of order T (where T is the width of the broader resonance when 
both resonances are separated) even in the general case of arbitrary number of open 
channels. This is clearly refiected in our numerical data presented in Fig. 1 b). 

Further progress can be made assuming that the quantum states localized 
in the chaotic regions of phase space obey random-matrix behavior. In this case, 
the width (W) becomes a statistical quantity and one can compute its probability 



distribution P(Ti) over a large number of anticrossings. Due to avoided crossings, 
this distribution shows a power-law decay P{Ti) ~ ^i^^'^ in the regime Fj <^ F, 
where F denotes the average width of the chaotic resonances. The derivation of the 
distribution -P(Fj) and a further discussion is defered to a forthcomming publication 
0. 

In summary, we have studied the lifetime and the phase space distribution 
of WG modes in a dielectric cavity. Dynamical tunneling leads to large fluctuations 
of the lifetimes as a function of the asymmetry parameter, as we demonstrated nu- 
merically for the open annular billiard. This system can be realized using dielectric 
microcavities. 
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